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Abstract. Recently, Ben Arous and Voiculescu considered taking the maximum of two free random 
variables and brought to light a deep analogy with the operation of taking the maximum of two inde- 
pendent random variables. We present here a new insight on this analogy: its concrete realization based 
on random matrices giving an interpolation between classical and free settings. 
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1. Introduction 





The free probability theory has been a very active field in mathematics over the last two decades, 
constructed in a deep analogy with classical probability theory. Nowadays, there is an unofficial dictionary 
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of concepts in both theories: many fundamental notions or results of classical probability theory, such as 
Law of Large Numbers, Central Limit Theorem, Gaussian distribution, convolution, cumulants, infinite 
divisibility have their precise counterpart in free probability theory. 

Recently, Ben Arous and Voiculescu [2] have added a new item in this dictionary, bringing to light the 
fact that the operation called the classical upper extremal convolution, which associates the distribution of 
the supremum X WY to the distributions of two independent random variables X and Y, has an analogue 
in free probability theory, called the free upper extremal convolution, which associates the distribution of 
the supremum X \/Y to the distributions of two free random variables X and Y. From the point of view 
of cumulative distribution function^, both of these operations have a concise interpretation: the classical 
(resp. free) upper extremal convolution of two probability measures with distribution functions F, G is 
the probability measure with distribution function FG (resp. max(0, F + G — 1)). The purpose of this 
paper is to construct random matrix models providing a better understanding of the free upper extremal 
convolution and a new insight on the relation between classical and free upper extremal convolutions. 

Let us first notice the following fact (which is an immediate consequence of remark |2.1[) : for /i, v 
probability measures on the real line, for (Xj)j>i and {Yi)i>\ two independent families of independent 
identically distributed random variables, with respective distributions fi, u, if for each N > 1, one de- 
notes by Zjv.i > Zm,2 > • • • > Zn,n the N largest elements of the multiset {X\, . . . , Xn, Y\, . . . , Y/v), 
then the empirical probability measure 53j=i fe N « converges almost surely to the free upper extremal 
convolution of fi and v as N tends to infinity. This result may appear as a coincidence. But, if we 
consider (Xi, . . . , Xn) and (Yi, . . . , Y/v) as the eigenvalues of two independent random matrices X( N > 
and YW which are invariant in law under conjugation by any unitary matrix, then (Zjv,i, ■ ■ ■ , Zn,n) 
are the eigenvalues of the supremurrQ of X^ and YW with respect to the spectral order introduced by 
Olson [5] and used by Ben Arous and Voiculescu in their paper. Since X W and F W are known to be 
asymptotically free, this gives actually a first interpretation of the free upper extremal convolution. 

We would like to give a deeper insight on the analogy between the classical and free settings. There is 
a morphism A v from the set of probability measures endowed with the classical upper convolution to the 
same set endowed with the free upper convolution. It maps any probability measure with distribution 
function F to the one with distribution function max(0, 1 + logF). Moreover, Ben Arous and Voiculescu 
have established that A v provides a remarkable correspondence between classical max-stable laws and 
free max-stable laws, which preserves the domains of attraction. In this paper, we shall give a concrete 
realization of this morphism via a random matrix model. 

More precisely, for each probability measure \i on the real line and each positive integer N, we shall 
define a law A n (/j.) on the set of N by N Hermitian matrices such that: 

a) a A N (^)-distributed random matrix is invariant under conjugation by any unitary matrix, 

b) for any pair /Lti,/X2 of probability measures on the real line and any pair M\,Mi of independent 
random matrices distributed respectively with respect to A N (\x{) and A N {jx-x), the law of the 
supremum M\ V M-i of M\ and is A]y(/Li), where fi is the classical upper extremal convolution 
of Hi and (12, 

c) if for all N, Mn is a random matrix with law A N (/i), then the empirical spectral law of Mm 
tends almost surely to A v (/x) as N tends to infinity. 

Such a model produces a new clue with regards to the relevance of the analogy between the two max- 
operations, and it may be used to provide more intuitive proofs of some of Ben Arous and Voiculescu's 
results. For instance, we have already mentioned that unitarily-invariant random matrices behave asymp- 
totically as free random variables. Therefore one can re-prove immediately, as a consequence of a), b) 
and c), the fact that A v is a morphism between the classical and the free upper extremal convolutions 
(and thus the formula Fxw = max(0, Fx + Fy — 1) for free random variables X, Y). We shall see later 

lr The cumulative distribution function of a random variable X with distribution fi is the function, denoted either by Fx 
or _F M , is defined by F^(x) = 00, x]), for all i£t 

2 The supremum of two Hermitian matrices is defined in section [2. jj I 
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that Ben Arous and Voiculescu's results about max-stable laws and domains of attraction can also be 
proved with our random matrix model. 

In a forthcoming paper, we will show how this approach is also appropriate to study the multivariate 
counterpart of the free upper extremal convolution. 

The paper is organized as follows. First, we shall focus on the upper extremal convolution for classical, 
free and Hermitian random variables. Then, we shall introduce our models, via a strategy which is close 
to the one of our previous papers [3] and [5], and set the main theorems. The last part will be dedicated 
to the proofs. 

2. The upper extremal convolutions 

2.1. For independent real random variables. Let X and Y be real independent random variables, 
with cumulative distribution functions Fx and Fy, and let us denote XVY = max(A, Y). The cumulative 
distribution function of X V Y is then equal to Fxyy = FxFy. 

2.2. For free random variables. Let (M, r) be a tracial Vy*-probability space, that is a von Neumann 
algebra M endowed with an ultraweakly continuous faithful trace-state r. The maximum X V Y of two 
self-adjoints elements X and Y in M is defined with respect to the spectral order introduced by Olson 
in 1971 (cf [9]; see also [HE]): 

• If p and q are two self-adjoint projectors, then q V q is the Hermitian projector on + 3(g). 

• If a and b are two self-adjoint elements, with resolutions of identity E a and then a V b is the self- 
adjoint clement with resolution of identity given by E a Vi?;,: a V b is the only self-adjoint element h such 
that for all real number t, 

X(t,oo)(h) = X(t,oo)(a) V X(t,oo)(b), 
where X(t,oo) denotes the indicator function of the interval (£, oo) and is applied via the functional calculus. 

In 2006, in [2J, Ben Arous and Voiculescu have determined the cumulative distribution function of 
XVY whenever X, Y are free: 

(1) V(Gl, FxvY{t)=max{0,F x (t)+FY(t)-l)=:Fx®Fy[t). 

It is easy to note that the function A v (u) = max(0, 1 + log it) is a kind of morphism between classical 
and free upper extremal convolutions: 

A V (F X F Y ) = A V (F X )UA V (F Y ). 

Notations: • In order to avoid the use of too many notations, we shall define the binary operation 
M on the set of probability measures on the real line with the same symbol as the operation U on the 
set of cumulative distribution functions. For all pair (j,, v of probability measures on the real line, /xM v 
is defined by: 

F,& V = F^UF U { = max(0, F„{t) + F v {t) - 1) ). 

• We shall also define the operator A v on the set of probability measures on the real line with the same 
symbol as the function A v (x) = max(0, 1 + log a;) on [0, +oo): for all probability measure [i on the real 
line, A v (/x) is the probability measure with cumulative distribution function: 

F A v M (t) = A v (F M (t)) ( = max(0,l + logF„(t)) ). 

Remark 2.1. The free upper extremal convolution also appears in a much simpler situation than the 
maximum of two free operators. Let X\, . . . , Xn, yi, ■ ■ ■ , UN £ and z\ > ... > z x the N largest 
elements of the multiset (a?i, . . . , xjy, yi, . . . , ?Mr). Then 53j=i * s the free upper extremal convolution 
°f jf Eili &xi and ± J2?=i S vi ■ 
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2.3. For independent Hermitian random matrices. Let N be a positive integer. The maximum 
operation introduced by Olson is also denned on the space of N x N Hermitian matrices, as a particular 
case of the definition given in section [2~2l 

Concretely, for A, B Hermitian matrices with the same size, A\/ B can be defined in the following way. 
Let {Ai > A2 > • • • > A p } be the union of the spectrums of A and of B. For i = 0, . . . ,p, let us define 
Ei = X]}=i(ker(A — Xjl) + kcr(_B — Xjl)). Then A V B is the Hermitian matrix with eigenvalues the X/s 
such that Ej ^ Ej—i, with associated eigenspace Ej fl (Ej_i). In other words, if, for any subspace F, Pp 
designs the orthogonal projector onto F, we have 

A V B = AiPej + A 2 P S2n ( Bj L) + hPE 3 n(Ei-) H ^ \> p e t c\(e£_ x )- 

The following proposition resumes the characteristics of the operation V on unitarily invariant random 
matrices. Recall that the empirical spectral law of a matrix is the uniform measure on the multiset of its 
eigenvalues counted by multiplicity. 

Proposition 2.2. If A and B are two N x N Hermitian random independent matrices, whose laws are 
invariant by conjugation by any unitary matrix, then 

• Ay B is an Hermitian random matrix, whose law is invariant under the conjugation by any 
unitary matrix; 

• the N eigenvalues of A\I B are the N largest of the 2N eigenvalues of A and B ( counted with 
multiplicity); 

• the cumulative distribution functions Fa, Fb and Favb of the respective empirical spectral laws 
of A, B and A\/ B are almost surely linked by the relation 

(2) F AyB = max(0, F A + F B - 1). 

Remark 2.3. It is noteworthy that formulas ([T]) and @ are identical, showing that the matricial and 
free upper extremal convolutions are somehow the same operation. As free random variables can be 
approximated by independent unitarily invariant random matrices, the rather trivial result @ provides 
an intuitive proof of formula ([!} . 

3. The matricial interpolation 

The aim of this section is to introduce random matrices giving quite a natural interpretation of the 
mapping A v . We simply follow the same approach as we did for the Bcrcovici-Pata bijection in [S] and 

Let Fp denote its cumulative distribution function for fi probability measure on the real line. For 

JL 

all integer k > 1, Ffi is the cumulative distribution function of a probability measure on the real line, 
which will be denoted by p,* . It is clear that if X\, . . . , are independent random variables distributed 
according to p^ , max{Xi, . . . , Xj-} should be distributed according to /j. 

Let us consider N, k > 1. Let X\, . . . , Xn be independent random variables distributed according to 
/ifc and U an unitary, Haar distributed, N x N random matrix independent of X\, . . . , Xjy. Define the 
random matrix M by 

(3) M = U 
Let Mi, . . . , Mk be independent replicas of M. 

Theorem 3.1. As k tends to infinity, the distribution of Mi V- • -VMfc converges to a probability measure 
A^f(p) on the space of N x N Hermitian matrices. This measure is invariant under the action of the 
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unitary group by conjugation and in the case where fj, has no atom, under A)y(/z), the joint distribution 
of the ranked eigenvalues is absolutely continuous with respect to /j,® n , with density 

at any point (ti, . . . , ijv) £ ^ N such that for all i = 1, . . . , N, F^iti) ^ (the density being set to zero 
anywhere else). 

Moreover, if ' v is another probability measure on the real line and A, B are independent random matrices 
respectively distributed according to A^(p,), A^(^), then A\/ B is distributed according to h v N (p), where 
p is the probability measure on the real line such that F p = F p F v . 

Corollary 3.2. With the preceding notations, we have 

(i) A^([i) is max- infinitely divisible, 

(ii) if n is max-stable, then so is A)(f(/x). 

In the following theorem, we are interested in the limit of the spectral measure of a h^ N (/^-distributed 
random matrix, when its dimension N tends to infinity. 

Theorem 3.3. Let fi be a probability measure on the real line. For all N > 1, let Mjv a random matrix 
with law A)y(/x), Ax, ■ ■ • , Ajv its eigenvalues and fiN — SiLi ^ its empirical spectral law. Then, when 
N goes to infinity, fiN converges weakly almost surely to A v (/i). 

The previous theorem may be used to provide new, more intuitive proofs of some of the Ben Arous 
and Voiculescu's results in [2]. 

A first example has been given in the introduction: it is a new derivation of the formula ([T]) defining 
the free upper extremal convolution in terms of cumulative distribution functions. 

For second example, the fact that A v maps any classical max-stable law to a free one is a direct 
consequence of (ii) of Corollary 13. 21 and of formulas (Q]) and @. 

Finally, let us consider the preservation of the domains of attraction. Ben Arous and Voiculescu have 
proved that for any classical max-stable law //, A v (/z) is freely max-stable and that for any cumulative 
distribution function F, for any sequences a^ > 0,bk £ K, we have 

(5) F(a k ■ +b k ) k — » F(a k ■ +b k ) ak — > A v (^(-)) 

k — >oo k — >oo 

(to state the reciprocal implication, one needs to choose a right inverse to A v , which we shall not do here). 
This result can be generalized as follows: for all sequence (F k ) of cumulative distribution functions, if F 
is a cumulative distribution function, we have 

(6) ^(-) fe — ► F(.) F k (-f k — > A v 0F (■)) 

k — >oc k — >oo 

(the proof of ([6]) is given in the last section of the paper). 

Let us now explain how our random matrix model gives an heuristic explanation of (JH). For sake of 
simplicity, let F k be equal to Fn , with F the cumulative distribution function of fi (for more general F k , 
a slight modification of our random matrix model is enough to apply what follows). 

For each N > 1 and for each k > 1, let us consider a family Mi(N, k), . . . , M k (N, k) of independent 
replicas of the random matrix M of (J3]), i.e. 

•.. 
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where X\ , . . . , Xjy are independent random variables with cumulative distribution function F <= and U 
an unitary, Haar distributed, N x N random matrix independent of X\, . . . ,Xj^. 



Let k tend to infinity. Following the definition of A^(fi) given in theorem 13. f I we get: 



(7) Mi(N, k) V ■ • • V M k (N, k) — ► M N in distribution, 

k — >oo 

where for all N, Mjv is a K y N (/i)-distributcd random matrix. 
Let now N tend to infinity. From theorem 13. 3[ we deduce: 

(8) Empirical Spectral Law(A/jv) — > A v (/x). 

AT— >oo 

• Now, let N tend to infinity before k does. According to the law of large numbers, the empirical 
spectral law of Mi(N, k) converges almost surely, for each k and i = 1, . . . , k, to the law with 
distribution function F*, which shall be denoted by /n* . Therefore, we deduce from the asymp- 
totic freeness of independent randomly rotated random matrices (more precisely by the equality 
between formulas fT]) and @), the following almost sure convergence: 

(9) Empirical Spectral Law(Mi(AT, ft) V ■•• V M k ( N, k)) — ► [i*® • • • M fi* . 

A— >oo 

J ommg 0, ©, ®, we get the following diagram 

k— >oo 

Mi (d, k) V ■ • • V M k (d, k) > (n) 

I I 

N -> oo N -» oo 

I I 
Empirical Spectral Law: Empirical Spectral Law: 

■ ■ -M fi* A V (/ i ) 
The right-hand side of ©, i.e. the fact that (^) ofc — ► A v (/i), only means that one can add an edge 

k — >oo 

— > between both bottom vertices of the diagram, i.e. that the operations k — > oo and N — > oo are 

— >oo 

commutative, which is quite expected. 

4. Proofs 

4.1. Preliminary results. Both of the theorems will be proved first for measures which arc absolutely 
continuous with respect to the Lebesgue measure, and then extended to all measures by approximation. 
In this context, the appropriate approximation tool is given by the following lemma. 

Lemma 4.1. Let us define, for ^ probability measure on the real line, the function F< _1> on (0, 1) by 
^f -1> («) = min{a; G K ; F M (x) > u}. 

(i) For \x,v probability measures on the real line, we have the following equalities (between quantities 
which can be infinite) 

inf{e > 0; F v (» - e) < < F„{* + e)} = \\F<- 1:> - F^W^ = mi\\X -Y^, 

where the infimum in the third term is taken on pairs {X, Y) of random variables defined on a same 
probability space with respective distributions fi, v. 

(ii) For any probability measure fi on the real line, for any e > 0, there exists a probability measure \i e 
on the real line such that F^ is smooth and 

F^C-£)<F^)<F^(. + e). 

Proof, (i) We prove these equalities by cyclic majorizations. 
- Let us prove first 

inf{e > 0; G R, F v (x - e) < F^x) < F v (x + e)} > \\F<- 1:> - F^'^W^. 
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Let us consider e > such that F v (» — e) < F (ti (») < F„(« + e). Since this is equivalent to F^(» — 
e) < F,(») < Ffj,(» + e), this inequation is symmetric in [i and v. Hence, it suffices to prove that 
F<- X> {u) - F<- X> {u) < e for all u e (0, 1). We have: 

F v {F<- x >{u)+e) > F,{F<- X >{u)) > u, 

hence, F<~ 1:> {u) < F<- 1;> (u) + e. 

- The inequality 

WF^ - F^Woo >ini\\X -YW^ 

is due to the fact that for U random variable with uniform distribution on (0, 1), X := F< ~ 1> (U), Y := 
F<~ 1> (?7) are respectively distributed according to \i, v. 

- To conclude, it suffices to prove 

inf \\X - riloo > inf{e > ; F,(. - e) < < F„(. + e)}. 

So let us consider a pair (X, Y) of random variables defined on the same probability space with respective 
distributions [i, v. Consider e > such that \X — Y\ < e uniformly on the probability space. Then, for 
all real number x, 

P(Y < x - e) < P(X < x) < P(Y < x + e). 

(ii) Let e > 0, and F a smooth non decreasing function on the real line such that for all k 6 Z, 
F(ke) = Ffj_(ke). Such a function exists obviously, as an example, one can consider F{x) = J_ f(t) dt, 
for 

fit) = - F (( k - - k£ ) 

fcez 

with smooth non negative function on the real line, with support contained in [—£,0] satisfying 
f_ E ip(t)dt = 1. Hence F is a cadlag non decreasing function on the real line such that F{x) tends 
to zero (rcsp. one) as x tends to — oo (resp. +oo). It follows that F — F^ E for a certain probability 
measure fi £ on the real line. Moreover, for all real number x, if k £ Z is such that ke < x < [k + l)e, 
then, one has: 

F^ix -e)< F^ c (ke) = F„(fce) < F^x) < F M ((fe + l)e) = F Me ((fc + l)e) < F^(a; + e). 

□ 

For any positive integer N, let us define max^r to be the function from U„>jvlR n to M. N which maps 
any vector x of K n , n > N, to the vector of the N largest coordinates of x ranked in decreasing order. 
The following property is a basic result, see for instance [7]. 

Proposition 4.2. Let F be the cumulative distribution function of a probability measure \i with no atom 
on the real line. Let, for n > N, X\, . . . , X n be independent random variables with law /i. The distribution 
of maxAr(Xi, . . . , X n ) has density 

(h, . . . ,t N ) i-> - N y h 1 >->t N F(t N ) n ~ 

with respect to n® N . 

Before stating the next proposition, we shall recall that for fi be a probability measure on the real line 
and k positive integer, is the probability measure on the real line with cumulative distribution function 

Ffi , i.e. the law m such that for X\, . . . , Xk independent random variables with law m, max(Xi, . . . , Xh) 
has law \x. 

Proposition 4.3. (i) Let fi be a probability measure on the real line. For all integer N > 1, the push 
forward of the probability measure (/jk)® kN nR kN by the function max converges weakly, as the integer 
k tends to infinity, to a probability measure on M. N denoted by /in. 
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(ii) When p has no atom, pm is absolutely continuous with respect to p® N , with density 

n i tl >-> t «ll^y 

at any point (t\, . . . , tjv) S such that for all i = 1, . . . , N, F^{ti) ^ 0. 

(Hi) Let us endow M. N with the norm \\x\\ = max^ |xj|. Then, for any probability measures p, v on the 
real line, 

inf -f 1 1 V — W^Hoo ; V,W random vectors defined on the same space with respective laws pn,vn} 
< inf{||X — y||oo i random variables defined on the same space with respective laws p, v\. 

(iv) For any pair p, v of probability measures on the real line, if p is the probability measure on the real 
line such that F p = F^F^, for all N > 1, p^ is the push-forward, by the function maxjv, of the probability 
measure pn ® vn on M. 2N . 

Proof. Let k,N > 1, and p a probability measure with no atom. From property 14. 2\ we infer first 

that pi = j:F^ 1 dp, and then that the distribution maxjv((/x^)® fcAr ) has a density with respect to p® N 
equal to 

(dk)\ , nW 4 



(UK)'. t— r 



at any point (ti,..., tjv) € such that for all i = 1, . . . , N, Ffa) ^ (the density can obviously be 
set to zero anywhere else). As k tends to infinity, this density stays uniformly bounded and converges 
pointwise to 

hence (i) for probability measures with no atom and (ii) are proved. 

Now, let us complete the proof of (i). Let us consider a probability measure p and, for all positive 
integer k, a family X(k, 1), . . . ,X(k, kN) of independent random variables with law p* . To prove (i), 
by Theorem 1.12.4 of [TU], it suffices to prove that for any real bounded Lipschitz function / on ~R N , 
the sequence E(/(maxjv(X(fc, 1), . . . , X(k, kN)))) converges as k tends to infinity. So let us fix such a 
function. We shall prove that the previous sequence is a Cauchy sequence. Let us fix e > 0. Let us 
consider p e as in (ii) of lemma 14. 11 Note that for all positive integers k, we also have: 

f|(.- £ )<f|(.)<F* (. + £). 

So, based on (i) of lemma [4~T1 we may suppose that for all k > 1, on the probability space where the 
X(k, z)'s arc defined, there is a family 

Y(k,l),...,Y(k,kN) 

of random variables such that 

JL 

(a) the Y(/c,i)'s are independent and distributed according to pi , 

(b) for all i = 1, . . . , kN, \X(k, i) - Y(k, i)\ < 2s almost surely. 

Note that p £ has no atom, hence by (a) and what we just proved, the sequence 

E(/(ma XA r(r(fc,l),...,y(fc,fciV)))) 

is a Cauchy sequence. Note also that by (b), if C is a Lipschitz constant for / with respect to the norm 
||x|| = maxi \xi\, then for all k, 

\E(f(m&x N (X(k, 1), . . . , X(k, kN)))) - E{f(mnx N (Y(k, 1), . . . , Y(k, kN))))\ < 2Ce. 
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Hence there is fco > 1 such that for all k,k' > ho, 

|E(/(ma Xjv (X(fc, 1), . . . , X(k, kN)))) - E(f(m&x N (X(k', 1), . . . , X(k', k'd))))\ < (AC + l)e. 

Thus we have proved (i) for any probability measure /i. 

Now, let us prove (hi). Consider fi, v probability measures on the real line. Note that by part (i) of 
lemma 14.11 it suffices to prove that for all positive e such that 

F v (»-e) <F M (.) <F„(« + e), 

for all a > e, there exists a pair V, W of random vectors defined on the same space with respective laws 
Hn,vn such that \V — W\ |oo < a. Let us consider such a positive e and a > e. For all k > 1, we also 
have: 

F$ (•-£)< F$ (•) < F$ (• + £)• 

So, according to part (i) of lemma |4~T| we shall consider, for all k > 1, a family 

X(k, 1), . . . , X(k, kN),Y(k, 1), . . . , Y(k, kN) 
of random variables defined on the same space such that 

(a) the X(k,i)'s arc independent and distributed according to , 

(b) the Y(k,i)'s are independent and distributed according to fi, 

(c) for all i = 1, . . . , kN, \X(k, i) — Y(k, i)\ < a almost surely. 

Let, for all k > 1, be the joint law, on M. 2N , of the random vector 

(max N (X(k, 1), . . . , X(k, kN)),max N (Y(k, 1), . . . , Y(k, kN)). 

The law of maxjv(X(fc, 1), . . . , X(k, kN)) (resp. of max^ (Y (k , 1), . . . , Y(k, kN))) converges weakly to /zat 
(resp. to i/jv) as k tends to infinity. It follows that the sequence is tight, i.e. relatively compact for the 
topology of weak convergence (Theorem 6.1 of [3]). Let r be the weak limit of a subsequence of Tfe. Let 
(V, W) be a r-distributcd random vector of M. N x M^. Then the law of V (resp. of W) is (resp. vn). 
Moreover, it is easy to notice that for all k, is supported by {(x,y) ; x £ M. N ,y € K , — y|| < a}. 
Indeed, let u,t 6 S„ such that X(k,<r(l)) > ■■■ > X{k,cr{kN)) and Y(k,r(l)) > ■•■ > Y(k,r(kN)); 
then for any i = 1, . . . , kN, 

X(k,<j(i)) = max mm(X(k, j). j G V) 

yc{i,...,fcw},#y= 4 v v ' JhJ 

> max min(Y"(fc, j) — a, j €E V) 

~ Vc{l,...,fcJV},#V=i v v y 

= F(fc, r(i)) — a a.s. 
By symmetry, this proves \X(k,<j(i)) — Y(k,r(i))\ < a a.s., and this implies 
\\max N (X(k,l),...,X(k,kN))-max N (Y(k,l),...,Y(k,kN))\\ = max \X(k, a{i)) - Y(k,T(i))\ 

i—l,...,N 

< a a.s. 

Letting k go to infinity, this establishes the inequality \\V — W\\ < a almost surely. 

To prove (iv), consider for any k > 1, 

X(k, 1), . . . , X(k, kN),Y(k, 1), . . . , Y(k, kN) 

is a family of independent random variables such that for all £ = 1, ... , kN, X(k, i) (resp. Y(k, i)) is 
distributed according to fi* (resp. v*). Then X(k, 1) V Y(k, 1), X[k, kN) V Y(k, kN) are i.i.d. with 
law yO*?, and consequently iin,vn an d are the respective weak limit distributions, as k tends to infinity, 
of 

max N [X(k, 1), . . . , X(k, kN)}, m&x N {Y(k, I),..., Y(k, kN)} 
(10) and max N [X(k,l)VY(k,l),...,X(k,kN)VY{k,kN)}. 
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Hence, by continuity and commutativity of the maximum operations, maxjv(/ijv (8> Vn) is the weak limit, 
as k tends to infinity, of the distribution of 

(11) max N [X(k, 1), . . . , X(k, kN),Y(k, 1), . . . , Y(k, kN)}. 

To conclude, we shall establish that the distributions of the vectors of (fTU|) and (fTT|) have the same weak 
limits as k tends to infinity. Obviously, it reduces to prove that the probability of the event 

{max N [X(k, 1) V Y(k, 1), . . . , X{k, kN) V Y(k, kN)} 

(12) ^ maxAr[X(/c, 1), . . . , X{k, kN), Y{k, 1), . . . , Y{k, kN)]} 

tends to zero as k tends to infinity. Note that this event is equivalent to the fact that there is i <G 
{1, . . . , kN} such that X(k, i) and Y(k, i) are both coordinates of 

max N [X(k, 1), . . . , X(k, kN), Y(k, 1), . . . , Y(k, kN)} 

and are not equal to any of the others X(k,j)'s and Y(k,l)'s. Hence, if one denotes by I{k) (resp. J(k)) 
the set of i's in {1, ... , kN} such that X{k, i) (resp. Y{k, i)) is one of the coordinates of 

max N [X(k,l),...,X(k,kN)} (resp. of max N {Y(k, 1), . . . , Y(k, kN)}), 

the event of (TT2")) is contained in the event I(k) n J(k) ^ 0, which probability is equal to , since we 

may choose for I{k) and J(k) to be independent and to have uniform distribution on the set of subsets 
of {1, . . . , kN} with cardinality N. Thus, (iv) is proved. □ 



4.2. Proof of proposition 12.21 The first statement is immediate. The third one follows immediately 
form the second one. To establish the second statement, it suffices to prove that almost surely, for any 
t £ R, Im(x( tjOCJ )(A)) + Im(x( t oo )(i?)) has the maximal dimension conditionally to the dimensions of 
Im(X(t,oo)(A) and Im (X(t,oo)(-B), i.e. that 

(13) dim[Im(x (ti00 )(A)) +Im(x(t,oo)(5))] = min{iV, dim[Im(x(t,oo)(-4))] + dim[Im(x( t ,oo)(-B)]} a.s. 

Let p = dim[Im(x(t,oo) (^4))]j € 1 = dim[Im(x(t,oo)(-B))]- Due to the unitarily invariance of A (resp. B), 
the law of Im(x(t,oo) (A)) (resp. Im(x(t,oo) (-£"))) is also invariant under the action of any unitary matrix, 
hence is uniform on the set of all subspaces of dimension p (resp. q). Moreover, Im(x({ j00 )(-<4-)) an d 
Im(x(t,oo) (B)) are independent. 

Let us consider a family (gij, i € {1, . . . , N} , j > 1) of independent complex standard Gaussian random 
variables. The matrix (gij, i, j G {1, . . . , N}) is distributed w.r.t. the Circular Unitary Ensemble. Hence, 
it is unitarily invariant and a.s. inversible. Therefore, the subspace generated by vectors {gij,i € 
{l,...,iV})f=i ( res P- (9%,j,i G {1, • ■ ■ ,N}) P jtl +1 , (gi,j,i e {l,---,^})^?) has dimension p (resp. q, 
min(iV,p + q)) a.s., and its law is uniform on the set of all subspaces of same dimension. This proves 
(HT 



Remark 4.4. Note that the proof of the second statement could also have been deduced from Theorem 
2.2 o/[6]. 

4.3. Proof of theorem [3TT1 Following proposition [2~2l if we denote = M x V • • • V M k , then 

is an Hcrmitian random matrix, whose law is invariant under conjugation with any unitary matrix and 

(k) (k) 

whose eigenvalues X\ > ... > are equal to the N largest of the Nk eigenvalues of Mi, . . . , Mk 
which are all independent with same distribution ^}/ k . Therefore, by theorem 4.3.5 of [8], there exists U 
unitary, Haar distributed and independent of the vector (A^ , . . . , A^ ), such that: 

( ^ 

M {k) =U 
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proposition 



. , A^ ) is the push forward of the probability measure 
kN by the function maxjv. Hence, the first part of the theorem follows from (i) and (ii) of 
and the second part follows from (iv) of the same proposition. 



Note that the distribution of the vector (A^ 



4.4. Proof of corollary 13.21 Both assertions are easy consequences of the last part of theorem 13. II 



(i) Let p > 1. By definition, F^ = [F 



Therefore A^(Fp) is the distribution of Hi V • • • V H p , 



with Hi, ... , H„ i.i.d. with law A^(F i ). This means that A)^(Fp) is max-infinitely divisible. 

(ii) Let p > 1. Since [i is max-stable, there exists a p > 0, b p G K such that fj, is the distribution of 
{flpXi + bp) V • • • V (dpXp + bp), where Xi, . . . , X p are i.i.d. with law fi. Let Hi, ... , H p be i.i.d. 
with law A)y(^). It is quite immediate from the construction that the image of the distribution 
of a p Xi + b p is the law of a p Hi + bblN where In denotes the identity matrix. Therefore, A)(r(/i) 
is the distribution of (a p Hi + bpl^) V ■ ■ • V (a p H p + bpljsi)- This proves that A^(fi) is max-stable. 



4.5. Proof of theorem 13.31 Suppose first that Fp is smooth on the real line. It implies that fi is 
absolutely continuous with respect to the Lebesgue measure and that d/i(x) = F^(x)dx. Since Fp is 
continuous on the real line, the almost sure weak convergence of the sequence of random probability 
measures {jln to A v (^t) is equivalent to the fact the for all rational number t, //jy ((— oo, i]) converges 
almost surely to max(0, 1 + logFp(t)). We shall prove that for all t S R, /t^r ((— oo,i]) converges almost 
surely to max(0, 1 + logF M (t)). Let us fix a real number t. 

The first step is to compute E [e^™^ -00 '*!)] : with the notations Aq = +oo, Xn+i = — oo, we have 



E 



\fi N ((-co,t}) 



= E 



N 



-P) 

e N l\ p+1 <t<\ r 



_p=o 



N 



EAlV-p 
e^^P(A p+1 <t < A p ) 

Now, note that for all real number T such that F(T) > 0, for all integer non negative m, 



OogfMt))' 



,d»(t) (logF^T))^ 1 



'(T,oo) ^mW 

Hence for p = 0, . . . , N — 1, 



m + 1 



and 



(-oo,T] 



Fp{t) m Mt) = - 



F M (T) m+1 



m + 1 



P(A p+ i <t < X p ) 



lei 



N N 1 



N 



t>t p+1 >--->t N H<t p <-"<t 



i=l ^ % ' P' 



And for p = N: 



ArjV p + OO i— logF M (i) JVI-jV 

'(* <A ^)= (]V3T)! i Fp{t N ) N (-\ogFp{t N )) N - 1 d i i{t N )= J -—— u ^e- N -du. 
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Thus, using the fact that J2p=o = &X $x°° (N~iy uN ~ le ~ U we § e * : 

s AAjv((-oo,t])l 



e 



p=0 p ' ^° 



(iV-1)! 



A f + oo 1 /— logF^tWjV 

J-Ne-V logF„(t) l iV _ -LJ 1 JO l iV _ 



'-e-^xog^w (Jv-i)i"" " ' Jo (N-iy/ 

The following inequalities will be useful, and they can be established from the Chernov inequality or 

proved directly: 

(14) 

y x e [0, 1), f J*" - u^e-^ du < x N e N ^~ x \ Vx > 1, [ + °° - u N ~ l e- Nu du < x N e N ^~ x l 

Jo \N — 1)\ J x (N — 1)\ 

We have three cases to consider: 

Case 1: -logF^t) > 1. If moreover A e [0, iVlog (- logF^t))), then -e _ & logjF^i) > 1, and by 



N N N-l e -Nu du 



e-TT log F„(i) (N - 1)! 
\+N(l-e—k ) log F M (t) / _x „ , A w JV(l+e"& logF„(i)] 

< e v / [— e ^ log F M (t) J e V / 

= e JV(l+log F„ (t)+log(- log F„(t))) < x 

since 1 + logF^t) + log(- logi^(i)) < 0. Hence for all N > 1, E [e A Aw((-oo,t])j < 2 , and from the 
Chernov inequality, we get then: for all e > 0, 

P{u N ((-00,t])> S} < inf 2e -eA = 2e -eAnog(-logF fl (t))) 

Ae[0,JVlog(- log F M (t))) 

This is enough, with Borel-Cantclli lemma, to prove the almost-sure convergence of oo, t]) to 

= max(0, 1 + log F(t)). 

Case 2: -log^(i) < 1. Then by [fTI]). 

— log F M (t) ^JV 

_ tt^-^-^du < (- log F^t)f e^'-W) = 
with <J t = 1 + logF M (t) + log(- log J^(*)) < 0. 
• Moreover, we suppose A > 0. Then 

,A[£jv((-oo,i])-(l+logF M (i))] 



e 



= e 

From the Chernov inequality, we infer then: for all e > 



< (e x F^(t) N e~ Ne ^ lo s F "(*) + e NS t j e -A(i+logF M (t)) 

iVlogF (i (t)(^-(l-e-^)) +eJ v( 5t -£(l+logF„(t))) 



>{A J v((-oo,i])-(l + logF,(i))> £ } < ^/ e -^(<^--*)) + ^M(i + v,^ 

^ e -Ar(logF^(t)(A-(l-e- A ))+A £ ) + e ^(5 t -A(l+logF M (t)+e))^ 



inf 

A>0 



2e- iVC 
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for some constant C > 0. This gives, using Borcl-Cantelli lemma: 
limsup oo, t]) < l + logF M (£) a.s. 

N >+oo 

• We suppose now A < 0. Using the same trick, we get 

'{/M(-0O,i]) - (1 + logF^)) < -e} < in f ( e -^(log^(t)(A-(l-e^))-Ae) +e iV(^-A(l+log^(t)- £ ))^ 

< 2e- NC ' 
for some constant C > 0. This gives: 



liminf /tjv((— oo, t]) > 1 + logF^i) a.s. 

N > + oo 



And the result follows: 



lim Aiv((— oo, t]) = 1 + logFfj,(t) — max(0, 1 + logF M (i)) a.s. 

W > + oo 

Case 3: logF AI (i) = — 1. Wc know that F^(u) tends to as u tends to — oo. Hence, there is u < t 
such that \ogF^{u) < — 1. For all N, /ijv(( — oo, i]) > /tjv((— oo, u]), hence, 

liminf /ijv((— oo, t]) > liminf /Ujv(( — oo, it]) = = 1 + logi^i) a.s. 

W > + oo AT >+oo 

Moreover, for any positive e, by continuity of F^ and by the fact that F^iu) tends to 1 as it tends to 
+oo, there is u > t such that 1 + log -F M (i) < 1 + log F M (it) < 1 + log F^(t) + e. Hence, 

limsup /ijv((— oo, t]) < limsup /ijv((— oo, it]) = 1 + logF M (it) < 1 + logF M (i) + e a.s. 

JV >+oo JV >+oo 

Thcrforc, lim^y > +oc /tjv((— oo, t]) = 1 + logF^i) almost surely. 

The theorem is proved in the case where F^ is smooth. Now, we consider a probability measure /x on 
the real line without making any assumption about F^. By Theorem 1.12.4 of [10j . the distance 



d(m, m ) := sup 



/" / dm — J f dm 1 



where the sup is taken on the set BL\ bounded Lipschitz functions / on the real line with Lipschitz 
constant < 1 and such that ||/||oo < 1, is a distance which defines the weak topology on the set of 
probability measures on the real line. 

Thus we have to prove that almost surely, as N tends to infinity, 

limd(/t w ,A v ( M )) =0, 

i.e. that for all e > 0, almost surely, for N large enough, 

(15) d(/t iV! A v ( / i)) <e. 

So let us fix e > 0. By part (ii) of lemma |4~T1 there exists a probability measure v on the real line such 
that F v is smooth and 

F v (* - e/6) < < F„(« + e/6). 

Note that the same obviously holds if one replaces [i by A v (/i) and v by A v (^), since for all probability 
measure m, F^ii m \ = max(0, 1 + logF m ). Hence, according to part (i) of lemma |4~T1 

(16) d(A v („),A v (i/))<e/6. 

Note also that by definition of A)y(/x), the vector (Xi > • • • > Xn) of ranked eigenvalues of Mjy has law 
Hn- Following part (hi) of proposition 14. 3[ one can suppose that for all N, on the same space as Mjy, 
there is a random matrix Nn with law A N (i/) with ranked eigenvalues (Yi > • • • > Y&) such that for all 
i = 1, . . . , N, \Xi — Yi\ < e/3 almost surely. If one denotes the spectral law of Nn by On, it implies that 

(17) VW > 1, d(/ijv, On) < e/3 almost surely. 
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Please note finally that by the first part of the proof, 

(18) almost surely, for N large enough, d(i>pf, A v (z^)) < e/3. 

Equations (fTB)) . (fT7|) . (JTHJ) , together, imply (fl"5j) . Thus, the theorem is proved. 

4.6. Proof of the implication ©. Let us end the paper with the proof of the implication ©. We 
consider a sequence Fk of cumulative distribution functions and a cumulative distribution function F such 
that at any x where F is continuous, Fk(x) k — ► F(x). Let us prove that for any such x, Fk(x) ak — ► 

k — >oo k — >oo 

A v (F(x)). Let us denote, for y S K, y + = max(0,y). Note that since for any a,b G R such that 6 < 0, 
we have (a + + = (a + by induction on k, we prove easily that for all G\ , . . . , distribution 
functions, 

GiM •••UGk = (Gi H \-Gk — k + 1)+. 

Thus we have to prove 

(1 + k(F k (x) - 1))+ — » (1 + logF(a;)) + , i.e. fc(F fc (.x) - 1) — > logF(a;), 

fc — >oo k — 'oo 

which follows directly from the hypothesis Fk(x) k — > F(x). 

k — ^oo 
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